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^ ■ Abstract 

m 

^ . We study cosmological solutions in the low-energy effective heterotic string theory, 

O, 

£ — which is the Einstein gravity with Gauss-Bonnet term and the dilaton. We show that 

O 

the field equations are cast into an autonomous system for flat internal and external 
^ . spaces, and derive all the fixed points in the system. We also examine the time evolution 

of the solutions and whether the solutions can give (transient) accelerated expansion 
of our four- dimensional space in the Einstein frame. 
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§1. Introduction 



The recent cosmological observations have confirmed the existence of the early inflation- 
ary epoch and the accelerated expansion of the present universe. ^ An important problem is 
then to derive such a model from fundamental theories of particle physics. The most promis- 
ing candidates for such theories are the ten-dimensional superstrings or eleven- dimensional 
M-theory, which are hoped to give models of accelerated expansion of the universe upon 
compactification to four dimensions. There are many attempts to derive such models, but 
most of them assume some additional matters or need special settings. From the viewpoint 
of the fundamental theories, however, it is desirable if such models are obtained without 
making special assumptions. 

It has been shown that a model with certain period of accelerated expansion can be 
obtained from the higher- dimensional vacuum Einstein equation if one assumes a time- 
dependent hyperbolic internal space 2 ) and that this class of models is obtained 3 ^ from what 
are known as S-branes 4 )' 5 ) in the limit of vanishing flux of three-form fields (see also Ref. 6)). 
For other attempts at inflation in the context of string theories, see, for instance, Refs. 7)-9). 
Unfortunately this class of models do not give sufficient inflation necessary to resolve the 
cosmological problems. 

On the other hand, it has been known that higher order corrections can give rise to 
inflationary solutions. 10 ) This is a very desirable setting since there are terms of higher 
orders in the curvature to the lowest effective supergravity action coming from superstrings 
or M-theory. n )~ 14 ) The simplest such correction is the Gauss-Bonnet (GB) term in the low- 
energy effective heterotic string. (We ignore other gauge fields and forms for simplicity.) It 
is thus important to examine what kind of time-dependent solutions are possible in these 
theories. 

There are many works discussing cosmology with the GB correction in four and higher 
dimensions (see, for instance, 15)— 19)). For example, it was shown that there are two expo- 
nentially expanding solutions in the higher-dimensional space, which may be called gener- 
alized de Sitter solutions since the size of the internal space also depends on time. 15 ) Note 
that this does not mean that the solutions gives accelerating expansion in four dimensions. 
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Another interesting claim is that it is possible to obtain an inflationary solutions if the coef- 
ficient of the Gauss-Bonnet term is negative, 17 -* which is not the case in the effective theory 
of the heterotic string, and hence may not be relevant in our consideration. Moreover most 
of the work considers pure GB term without dilaton or assumes constant dilaton, which is 
not the effective theory of the heterotic string, and does not discuss cosmological solutions 
with dynamical dilaton in higher dimensions. It is thus important to analyze the system 
including the dynamical dilatons. Some attempt to obtain inflationary solutions in M theory 
with higher order quantum corrections has also been made. 20 ) 

Recently a more interesting approach is considered for Einstein theory with some ad- 
ditional scalars. 21 )' 22 ) In this dynamical system method, one considers the solution space 
restricted by the constraint equation resulting from a component of the Einstein equation. 
If the field equations are written as an autonomous system, we can find fixed points in this 
space. Then all possible solutions are expressed as trajectories between these fixed points 
in the solution space. This is a very powerful method to examine possible solutions which 
is applicable even if exact solutions are not available. In particular, it is possible to find 
solutions with (transient) accelerating expansion which may be relevant to cosmology. In 
fact, the existence of an eternally accelerating solution, first found in Ref. 8), is established 
for hyperbolic internal and external spaces without giving explicit solution. 21 ) 

In this paper, we consider cosmological solutions with a dilaton field and the GB cor- 
rection from heterotic string theory by extending the above dynamical system method. We 
find that the field equations may be cast into an autonomous system for flat internal and 
external spaces for both theories with and without dynamical dilaton. We derive all the fixed 
points and analyze their stability in the system. We also examine the time evolution of the 
solutions and investigate whether the solutions can give (transient) accelerated expansion of 
our four-dimensional space in the Einstein frame. 

This paper is organized as follows. In § 2, we first write down the action of the Einstein 
and GB theory, and our metric for D-dimensional space. We then summarize the field 
equations. In § 3, we analyze the theory without the dilaton and find the solution space and 
accelerating solutions. We show that the field equations become an autonomous system, and 
find fixed points. We can see how the solutions evolve in time by looking at the solution 
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space and the flow. We find that there is one fixed point corresponding to expanding solution 
with acceleration, but it gives a singular super-inflation. In § 4, we extend the analysis to 
the theory with the dilaton. At first sight, the field equations cannot be reduced to an 
autonomous system, but judicious choice of the time variable enables us to do it. We 
then discuss their flow and properties of the fixed points of our system. § 5 is devoted to 
conclusions. 

§2. Field equations 



We consider the low-energy effective action for the heterotic string: 

1 



2k? d 



d D x^ 



R + A(d^Y + d 2 R 2 GB \ , 



(2-1) 



where /^-dimensional gravitational constant, is a dilaton field, a 2 = a'/ 8 is a 

numerical coefficient given in terms of the Regge slope parameter, and R GB = R llV p a R tll ' pc ' — 
AR^R^ + R 2 is the GB correction. In the Einstein frame the dilaton is minimally coupled 



to the metric and has a canonical kinetic term 
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2k 2 d 



d D Xy 



R 



a 2 e-^R 2 GB 



where 9fiu = e" 4 ^- 2 ^, = ^/8/{D 
metric in D-dimensional space 



2) and 7 = y/2/(D-2) 



ds 2 D = -e 2u ^dt 2 + e 2u ^ds 2 p + e 2u ^ds 2 q , 



(2-2) 

Let us consider the 
(2-3) 



where D — 1 + p + q. The external p- and internal g-dimensional spaces (ds 2 and ds 2 ) are 
chosen to be maximally symmetric with the signature of the curvature given by a p and a q , 
respectively. Though we are mainly concerned with flat internal and external spaces in this 
paper, it may be useful to give field equations for more general case. 
We find that the Riemann tensors are given by 

R t itj = e 2u °Xgij , 
R'at^e-^Yg^, 



R 



jki 



' 2uo A p (g\ gjl - g\g jk ) , 
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R\jb = e 2u °uiU2g l jg a b , 

R\cd = e- 2uo A q (g a c g bd - g a d g bc ) , (24) 
where i,j and a, b run over p- and g-dimensional spaces, respectively, and 

A p = til + v P e 2(u °- ul \ A q = ii 2 2 + a q e 2(u °- U2 \ 

X = Hi — UqUi +ii 2 , Y = u 2 — UqU 2 + u 2 , . (2-5) 

The GB term is given by 

Rqb = e" 4M0 {ps^5 + 2p 1 q 1 A p A q + q 3 A 2 q + Ainu 2 {p 2 qA p + pq 2 A q ) + Ap x q x u\u 2 2 
+ ApX [{p - 1) 2 A P + qi A q + 2(p - l)gui« 2 ] 

+ 4gF [p^ + (g - 1) 2 A ? + 2p(q - l)inu 2 ] }, (2-6) 

where we have defined 

(p — m) n = (p — m){p — m — l)(p — m — 2) • • • (p — n) , 

(g - m) n = (g - ra)(g - m - l)(q - m - 2) • • • (q - n) , (2-7) 

Multiplying (2-6) by \f~-geT 1 ^ = e U0+pui+qU2 ~' y ^ and making partial integration, one finds 
that the action reduces to the following (up to an overall factor): 
(1) Einstein-Hilbert action 



£ _ e ~u +pu 1 +qu 2 

(2) GB action 



PiA p + q x A q - 2(pi«i 2 + pqu x U2 + qiu 2 2 ) + 



£ 2 = e -3u 0+P u 1+q u 2 -y^ p3Ap 2 + 2piqi A q A p + q 3 A q 2 

- AA p (pzU 2 + p 2 qu x u 2 + piqiii 2 2 ) - AA^p^u 2 + pq 2 u x u 2 + g 3 ^ 2 2 ) 
+ ^(2p 3 Mi 4 + 2p 2 qii 1 3 ii 2 + 3pigi«i 2 w 2 2 + 2pq 2 u x u 2 + 2q 3 u 2 ) 



+ 470 



(P2«i + Piqu2)A p + (pgi«i + q2U2)A q - (p^? + g 2 « 2 )] | . (2-9) 



If we set = 0, this agrees with the results in Ref. 20). 



Now the field equations are 

F = F 1 + F 2 = 0, (2-10) 
F (P) = fir) + f (p) + x + + y ^(p) + h (P)^ _ z t ( P ) = o j (2 . n) 

= A ( 9 ) + /a ( 9 ) + y + ^) j + x + ftfa)) _ Z ita) = o , (2-12) 

F = Z + a 27 e 2u °-^ B = , (2-13) 

where Rq B is given in Eq. (2-6) and 

Z = (f) + (—u + piii + qu 2 )^> , 
Fi =p 1 A p + q x A q + 2pqii 1 u 2 - ^0 2 , 

= (p - 1) 2 A P + + 2(p - 1)^ 2 + , 

= p^ + (q- l) 2 A q + 2p(g - 1)^ 2 + ^0 2 , 

^ = 2(p - 1) , 

9 { q) = 2(q - 1) , 

/4 P) = 2g, 

h ( f=2p, (2-14) 



and 



F 2 = a 2 e- 2u °-^|p3-4 2 + 2p 1?1 A p 4 z + q 3 A 2 q + A(p 2 qA p + pq 2 A q + p^inu^u^ 
- 470 [(p 2 ui + piqu 2 )A p + (pqiiii + q 2 ii 2 )A q + 2(p 1 qii 1 + pgiw 2 )wi« 2 ] | , 

/W = a 2 e- 2u °-^{(p-l) A Al + 2(p-l) 2qi A p A q + q 3 A 2 q 

+ 4 [(p - l) 3 gA p + (p - l)g 2 A 9 + (p - l) 2 giMiw 2 ] wi« 2 
+ 4 7 0[((p - 1) 2 A P + + 2(p - l)gui«2)(ui + 70) 
+ 2((p - l) 2 UiA p + q x u 2 A q + (p - l)guitt2(ui + u 2 ))] j, 

/<«> = a 2 e- 2 "»-^{p 3 ^ + 2 Pl ( g -l) 2 ^ + ( g -l) 4 ^ 

+ 4 [p 2 (g - + p(q - l) 3 A q +pi(q- l) 2 UiU 2 ] U\U 2 
+ 4^[( Pl A p + (q- l) 2 A q + 2p(g - l)uiu 2 )(u 2 + 7 0) 
+ 2(piMiA p + (g - 1) 2 m 2 A (? + p(g - l)uiii 2 (ui + w 2 ))] | , 
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4(p - l)a 2 e 



(p - 2) 3 A P + giA, + 2(p - 2)qii 1 u 2 - 2 7 ((p - 2)in + gw 2 )c/> 



4(g - l)a 2 e- 2u °-^ Pl A p + (q - 2) 3 A q + 2p(q - 2)uiti 2 - 2 7 (pui + (g - 2)u 2 )0 



4ga 2 e 



-2uo-70 



(p - l) 2 A p + (g - 1) 2 A ? + 2(p - l)(g - l)uiu 2 



- 2 7 ((p - l)«i + (g - l)u 2 )^ 



= 4pa 2 e 



-2m -7<j!) 



;(p) 



juv — 



(p - l) 2 A p + (q- l) 2 A q + 2(p - l)(g - l)uiu 2 
-27((p-l)ui + (g-l)« 2 )^], 
a 2 e- 2u °-^4 7 [(p - l) 2 A p + Ql A q + 2(p - l)g^ 2 , 



;(<z) 



= a 2 e- 2 "°-^4 7 Pl A p + (q - l) 2 A q + 2p(g - 1)^ 2 



(2-15) 



The basic Eqs. (2-10) - (2-13) are not all independent. They satisfy 
F + (pui + gw 2 - 2m )F = puiF^ + qii 2 F^ - 0F . 



(2-16) 



We are now going to examine cosmological solutions in this system. In this paper, we only 
consider flat internal and external spaces, i.e., a p = a q = 0. Henceforth, we set p = 3 and 
q = 6 though we write formulae for more general cases as much as possible. 

§3. Solutions in Einstein and Gauss-Bonnet theory 

In this section, let us first consider the theory without dilaton (or the case when dilaton 
is constant) as a consistency check. Namely we set <fi = by hand to investigate the system 
without the dilaton field, and study possible cosmological solutions in Einstein and Gauss- 
Bonnet theory without dilaton. This is the system examined in Ref. 15), but the question 
whether the accelerating expansion occurs or not in the four-dimensional spacetime was not 
examined, and we clarify this point by making systematic analysis of the solutions by the 
dynamical system method. We can set uq = by using time-reparametrization invariance. 
It is also possible to put a 2 = 1 by choosing a suitable unit of time. 20 -* 

Equation (2-10) is a constraint equation, and any cosmological solutions must satisfy 
this. In this sense, this gives the space in which all the possible solutions live, which we call 
solution space. This is depicted in Fig. 1 in (ui, ■u 2 )-plane. 

Solving Eqs. (2-11) and (2-12) for Hi and u 2) we find that the field equations (2-11) and 



7 




-4 -2 




Fig. 1. Solution space and flow. The dots indicate fixed points. 



(2T2) become an autonomous system for u\ and ii2- We then find the five fixed points of 
these variables for p = 3, q = 6 in the unit of a 2 = 1 are given by 



which are also shown in Fig. 1. We can also derive the flow of the solutions along the time 
lapse between the fixed points as shown in the figure. Due to the time-reversal symmetry of 
the system, the figure is symmetric under ir rotation (with the reversed time flow). All this 
agrees with the results in Ref. 15). 

Our cosmological model is higher-dimensional, and there are two kinds of frames that we 
can take to discuss cosmologies, the original frame and the Einstein frame in four dimensions. 
Note that this Einstein frame is different from the one defined in Eq. (2-2) with respect to the 
dilaton. Instead it is a new frame which is defined to eliminate the scalar fields which appear 
from the internal space by Kaluza-Klein compactification to the external space. This is the 
frame in which the Newton constant is really constant. We must determine which frame is 
important for a successful inflationary scenario. Since flatness and horizon problems should 
be explained in our four-dimensional spacetime, that is, in the Einstein frame, we should 
require a successful inflation in the Einstein frame. 

Now let us examine if there is any region where the accelerating expansion is realised in 
the four-dimensional Einstein frame. The Einstein frame is obtained by 




(0,0), (±0.88603, T-0.13845), (±0.48296,^0.34141), 





2q 

= e p- 1 



ds% + e 2u2 ds 2 q . 



(3-2) 
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So 



ds 2 E = e^ U2 (-dt 2 + e 2ui ds 2 p ) 



-dr 2 + a 2 (r)ds 2 (3-3) 



where we have defined the cosmic time r and scale factor by 



dt 

For p = 3, q = 6, the condition for expansion is 

da dt da 



^L = e ^ i a(r ) =e «i+^r« 2 . (3.4) 



, , ( (in + 3« 2 )e ui > 0, (3-5) 
dr dr dt 



and the condition for accelerated expansion is 

d 2 a dt d 



= {iix + 3w 2 + («! + 3u 2 )ui}e Ul - 3u ' 2 > 0. (3-6) 

Substituting u\ and «2 into (3-6), we find that the accelerating regions are those shown in 
Fig. 2, and the solution space in these regions are depicted in Fig. 3. 

We have also examined the stability of the fixed points. We find that the fixed points, 
(0, 0), corresponding to the flat Minkowski space, (-0.88603, 0.13845) and (0.48296, -0.34141), 
corresponding to a contracting universe, are unstable. The fixed points (—0.48296, 0.34141) 
and (0.88603, —0.13845) are stable. Among these, only the last one gives accelerating ex- 
pansion. The behavior of the scale factor is like a oc |r| -1 - 13 for negative r. This is what is 
called super-inflation and exhibits singularity near r ~ 0. There is a solution flowing into 
this fixed point which exhibits accelerating expansion for its whole evolution. Such a super- 
inflation may also arise in phantom cosmological models. 23 ^ This is also called a Big Rip 
and should be avoided. 20 ^ However, there are several (transient) accelerating cosmological 
solutions. For example, there is a solution coming out of the flat Minkowski space flowing 
into the direction of positive u\ and negative iii with transient acceleration. Whether this 
solution gives viable cosmological solution or not remains to be examined. 
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Fig. 2. accelerating region in the solution space. 
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Fig. 3. accelerating solutions. 



§4. Solutions with a dynamical dilaton 

In this section, we extend our analysis to the more interesting theory with a dynamical 
dilaton with p = 3 and q = 6, which appears as a low-energy effective theory of the heterotic 
string. We note again that uq = can be chosen by time reparametrization, and the choice 
of a suitable unit of time can be used to set a.2 = l- 20 -* 

In this system, there are exponential factors of the dilaton in the field equations (2-10) - 
(2-13), and this appears to prevent us from writing them as an autonomous system. However, 
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if we introduce new time variable T by 

d t = e*'*a r , i.e. ^ = e^ 4 , (4-1) 

then it is possible to rewrite them as an autonomous system. In what follows, derivatives 
with respect to T will be denoted by the prime '. Then the field equations (2-10) - (2-13) 
remain the same if we make the following replacement in Eqs. (2-5), (2-14) and (2-15): 

iii^u'l + ^uW, u 2 ^ u 'l + ±u' 2 <f>', 0^0" + I(0') 2 , (4-2) 

and remove the exponential factors in Eqs. (2-13) and (2-15). 

This is again an autonomous system for x = u[, y = u' 2 and z = 0'. Among these, 
the constraint (2-10) gives the solution space. In this case, because we have 3 variables 
x, y and z, the space consists of 2-dimensional surfaces embedded in 3 dimensions. The 
surfaces have the shape of hyperbolic surfaces. Since it does not seem to be so instructive 
to draw the surface in 3 dimensions, we show the shapes of slices of the solution space at 
4>' = 2,1,0.7,0.585906,0.3,0 in Figs. 4-9, respectively. We see that reconnections of the 
surfaces occur as 4>' varies. Note also that the region for <fi' < has just the 7r-rotated shape 
due to time reversal symmetry. 

We find that there are seven fixed points in this system 

(x, y, z) = M(0, 0, 0), Pi(^0.292373, ±0.36066, ±0.954846), 

P 2 (±0.91822, ^0.080285, ±0.585906), P 3 (±0. 161307, ±0.161307, ^9.30437), (4-3) 

where the labels are indicated for upper signs and those lower signs are denoted with tildes. 
Their properties are summarized in Table I. 

The expansion criterion (3-5) with replacement (4-2) tells us that the solutions of Pi, 
P 2 and P 3 give the expanding solutions in the Einstein frame. Among these, only P 2 gives 
accelerating expansion. In this accelerated solution, we have 

T - —pi 1 — - p 3«2-4>/4 _ -0.387T (A a\ 

$ ' ' dT ■ 1 } 

We then find that this solution with 

a(r) = e Ul+3M2 = e a677T ~ |r|- L75 , (4-5) 
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Fig. 6. Slice at <f>' = 0.7 in the solution space. Fi S- 9 - Slice at <f>' = in the solution space. 

gives again a super- inflation and r changes from -oo to as T changes from — oo to oo. 

In order to study the stability of the fixed points, we substitute linear perturbations 
x — > x + 5x, y — > y + Sy and z — > z + 5z about the fixed points into the field equations (2-10) 
- (2-13). To the first order in the perturbations, we obtain two independent equations of 
motion which can be written as 




where M. is a 2 x 2 matrix. Stability requires that both the eigenvalues of the matrix A4, 
Ai and A2 be negative. Our analysis shows that M, Pi, P2 and P 3 are unstable while Pi, P2 
and P3 are stable. 
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Label 


(x,y,z) 


Eigenvalues (Ai, A 2 ) 


Stability 


da 
dr 


d 2 a 
dr 2 


M 


(0,0,0) 


(0,0) 


unstable 


- 


- 


Pi 


(0.292373, -0.36066, -0.954846) 


(1.52555,1.52555) 


unstable 


< 


< 


Pi 


(-0.292373, 0.36066, 0.954846) 


(-1.52555,-1.52555) 


stable 


> 


< 


P 2 


(0.91822, -0.080285, 0.585906) 


(-2.41943,-2.41943) 


stable 


> 


> 


P 2 


(-0.91822, 0.080285, -0.585906) 


(2.41943,2.41943) 


unstable 


< 


> 


P 3 


(0.161307, 0.161307, -9.30437) 


(0.874329,0.874324) 


unstable 


> 


< 


P 3 


(-0.161307, -0.161307, 9.30437) 


(-0.87433, -0.874324) 


stable 


< 


< 



Table I. Fixed points of the autonomous system and their properties. 



The flow diagram for solutions around the fixed points is drawn in Fig. 10. We can 
use it to examine what kind of solutions are possible. For example, there are solutions 
starting from a decelerated expanding region which approach the accelerated expanding 
solution (P2), solutions starting from a decelerated contracting region which approach the 
accelerated contracting solution (P3), and solutions starting from a accelerated expanding 
region which approach the decelerated expanding solution (Pi). We see from this figure that 
there are several accelerating cosmological solutions in this theory including those flowing 
into non-accelerating fixed point Pi and those flowing into P2 with Big Rip singularity. It 
is possible that stringy effects resolve this kind of singularity and these solutions may give 
viable cosmologies. 

It is interesting to investigate whether these solutions give viable cosmological solution or 
not. A step towards this is to examine if we can get enough e- folding for solving cosmological 
problems. A preliminary investigation of the solution flowing into the fixed point P2 indicates 
that it is hard to get enough e-folding number before arriving at the fixed point but we can 
easily get sufficient e-folding number if the solution arrives at the fixed point. 

When we consider the accelerating expansion of the present universe, the fine-tuning 
problem is always a nagging problem. To partially answer this question, we have examined 
solutions by changing initial conditions near the fixed point P 2 , and find that there are 
several solutions flowing into P 2 , as shown in Fig. 10. This means that there are certain 
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Fig. 10. Solution space and flow in the case with a dynamical dilaton. The solid (red) lines 
correspond to cPa/dr 2 > and the dashed (green) lines correspond to d 2 a/dr 2 < 0. 

range of initial conditions which lead to the accelerating expansion. In this sense, these 
solutions have the possibility of explaining naturalness of the accelerating expansion. To 
examine how large area of these initial conditions can give such a behavior and whether the 
present model can give realistic one need further study. 

§5. Conclusions 

In this paper we have investigated cosmological solutions in the Einstein theory with 
GB correction with and without a dilaton in higher dimensions. We are interested in this 
theory because this is the low-energy effective theory of the heterotic string, and examined 
what solutions are possible by the dynamical system method. For flat internal and external 
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spaces, we have shown that the field equations can be written as an autonomous system for 
both the theories with and without dynamical dilaton. We obtained the fixed points and 
analyzed their stability. We have found that both in the GB correction with and without 
dilaton, there are solutions with accelerating expansion. Some of them are super-inflation 
with future singularity. 

The analysis in Ref. 21) indicates that even if there is no interesting cosmological solution 
in the Einstein theory for flat internal and external spaces, there may exist an interesting 
solution with eternally accelerating expansion for curved spaces. The existence of such 
a solution was originally suggested in Ref. 8) by a perturbation around non-inflationary 
solution, and it was shown that the solution is eternally expanding with acceleration after 
some time. However, due to the limitation of the perturbation, the detailed properties of the 
solution (like its eternal accelerating property for whole time) was not clear. The powerful 
method of dynamical system allowed to show that the solution is eternally expanding with 
acceleration for the whole time. 21 ) It is thus possible that similarly interesting solutions may 
exist in our Einstein-Gauss-Bonnet gravity coupled to dilaton. It would be very interesting 
to extend our analysis to curved external and internal spaces, and check if there may be 
additional interesting solutions. 
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